A solvable model of fracture with power-law distribution of fragment sizes 
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The present letter describes a stochastic model of fracture, whose fragment size distribution can 
be calculated analytically as a power-law-like distribution. The model is basically cascade fracture, 
but incorporates the effect that each fragment in each stage of cascade ceases fracture with a certain 
probability. When the probability is constant, the exponent of the power-law cumulative distribution 
lies between —1 and 0, depending not only on the probability but the distribution of fracture points. 
Whereas, when the probability depends on the size of a fragment, the exponent is less than — 1, 
irrespective of the distribution of fracture points. 

PACS numbers: 02.50.-r, 46.50.+a, 05.40.-a 



Fracture and fragmentation are ubiquitous in nature, 
and the comprehension and control of fracture are very 
important not only in science and engineering but in 
our daily lives. Numerous experiments have confirmed 
that fragment size distributions mainly follow power-law 
distributions [lHS], but other types, including lognormal 
distributions (J-Q, have been also observed. From a the- 
oretical point of view, various models have attempted 
to derive power laws pTT-flOj . but there seems to be no 
decisive model which briefly and analytically explains a 
power-law distribution of fragment sizes without using 
specific breaking mechanisms. Some stochastic models 
[Ill4l3| | can successfully provide power-law distributions, 
but they cannot be applied directly to fracture phenom- 
ena. 

A lognormal distribution of fragment sizes can be ex- 
plained by a quite simple model of cascade fracture [lij . 
In this model, one rod of length L breaks into two frag- 
ments at a randomly chosen point, and each of the two 
fragments again breaks into two sub- fragments, and so on 
(see Fig. [I] (a)). The length of one of the fragments af- 
ter the n-th stage of fracture is expressed as £1^2 ■ • ■ £rA 
where £i,£2i " " " i£n are random numbers between and 
1. This process is referred to as 'multiplicative', because 
the length of a fragment is given by multiplying the pre- 
vious length by £j. Assuming that £1, • • • , £ n are indepen- 
dently and identically distributed, and that the variance 
of log£i is finite, one can prove by the central limit theo- 
rem for log £j that the fragment size distribution exhibits 
a lognormal distribution when b>1. 

In the present letter, we slightly modify the above mul- 
tiplicative model of cascade fracture so that the resulting 
fragment size distribution becomes power-law like. 

Our model also starts with one rod of length L, and 
fragments repeatedly break into two sub-fragments. A 
fracture point is given by a random number £ G (0, 1) 
drawn from a probability density function <?(£). The 
difference from the above simple multiplicative model is 
that each fragment ceases fracture with a constant prob- 
ability p, which we call the "stopping probability" (see 
Fig. Q] (b)). Whether each fragment stops fracture or 
not is determined independently; once a fragment ceases 
fracture, it never restarts fracture any more, and we call 



such a fragment "inactive" . 

In order to analyze this model, we focus on the cumu- 
lative number Nr,{x) of fragments, which represents the 
expected number of fragments larger than x. The initial 
length L of the rod is specified as a parameter. Ni,(x) 
satisfies the following equation. 

N L {x) = (1 - p) [ {N £L (x) + N (1 _ 0L (x)} <?(£)d£ + p. 



(1) 

The first term (integration) at the right hand side repre- 
sents that the initial rod breaks into two fragments of the 
lengths £L and (1 — £)L with probability (1 — p), and the 
last term represents that the initial rod becomes inactive 
with probability p. 

If we rescale our length scale by a factor a(> 0) 
and observe fracture processes, the length of the initial 
rod is aL in the new scale, and the cumulative number 
Nl(x) corresponds to N a r J (ax). This seeming difference 
is just brought by rescaling, so we have a scaling relation 
N a L(ax) — Nl(x) or N a L(x) = Nl(x/o>). Using this 
relation to convert all subscripts in Eq. ([1} into L, and 
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FIG. 1: (a) The simple model of cascade fracture, where the 
fragments at each stage break into two sub-fragments. The 
resulting fragment size distribution is a lognormal one. (b) 
Our proposed model, which is similar to (a) except that each 
fragment ceases fracture with a constant probability p. 
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introducing 



N L (x) = N L (x) + 



1 - 2p 



(2) 



in order to eliminate the last term "p" in Eq. ((T|), we 
obtain a homogeneous equation of Nl' 



N L (x) = (l-p) j {N L 



Nt 



We assume a power-law form Nl(x) — Cx °, where C 
and /3(> 0) are both independent of x. Then, we have 



/ {^ + (l-0*}ff(0d£ 



1. 



(3) 



The exponent 8 is determined by this equation; hence 8 
generally depends on both p and g. Furthermore, it is 
natural to assume that the fracture is left-right symmet- 
rical if a rod is uniform, that is, g(£) = g(l — £) for any 
£ G (0,1). This condition simplifies Eq. (|3|) to 



2(1 -p) I ^9(0^ = 1 



(4) 



We finally determine the coefficient C. By the defi- 
nition of the cumulative number, Nl(L) is the expected 
number of fragments larger than L, and it is equal to 
the probability with which the initial rod ceases fracture. 
Thus, N L {L) = p. On the other hand, N L (L) = CL^ 13 
is immediately obtained. Therefore, it follows by Eq. 
that 



C = 



2p(l - p) 



1 - 2/5 

Eventually, the complete solution is 



N L (x) = 



1 -2p 



l-2p 



M* 1 (!)"'-'}• 



coupled with Eq. ([3]) for the determination of f3. Note 
that Nl is an exact power law, but Nl is not exactly be- 
cause of the presence of the second term "—1". Nonethe- 
less, Nl(x) can be approximated by a power law if the 
second term is negligible, i.e., 2(1 — p) 3> (x/L)P, or if x 
is sufficiently smaller than L and p is also small. 

The solution /3 of Eq. ^ or Q cannot be expressed 
explicitly for general probability density g. We provide 
two examples of calculations of 8. The mathematically 
simplest instance is g(£) = 5(£ — 1/2), where S is the 
Dirac delta function. In other words, the fracture points 
are at the middle of the fragments. Equations (|3J) and 
Q in this case are both reduced to 



2(1 



= 1, 



and the solution is j3 = 1 + log(l — p)J log 2. 

In the second example, the fracture point is distributed 
uniformly over each fragment, i.e., g (£) = 1 for all £ G 
(0, 1). Equation Q becomes 



2(1 -p) 



/3 



= 1, 



and the solution is /? = 1 — 2p. 

Both in these two specific examples, (3 = 1 + log(l — 
p)/log2 and /3 = 1 — 2p, the possible ranges of (3 and p 
are restricted by each other: p < 1/2 follows from /3 > 0, 
and (3 < 1 follows from p > 0. Namely, the reasonable 
value of the stopping probability is0<p<l/2 and the 
possible value of f3 is within < j3 < 1. 

The same restrictions for p and j3 also hold for a gen- 
eral probability density g. For a fixed g, we consider 
f(fi) ■= Jq 1 + (1 - £)£} g(0dC, which is defined at 
least in f3 > 0, and continuous and differentiable. Clearly, 
Eq. §S§ is expressed as f((3) = 1/(1 — p). Using the 
normalization Jq g(0^ = 1j we bave /(0) = 2 and 
/(l) = 1 for any g. Hence, the intermediate value the- 
orem in elementary calculus insures that the equation 
/(/?) = 1/(1 — p) has at least one solution j3 G (0, 1) if 
1 < 1/(1 - p) < 2 (or < p < 1/2). On the other hand, 
/(/?) is decreasing because 



f(P) 



/V iog£ + (i - io g (i - e)}5(C)de < o, 

Jo 



therefore the correspondence between 8 and p is one- 
to-one, i.e., 8 can be determined uniquely for given p. 
Positivity p > holds in /3 < 1, and 8 > holds in 
1/(1 - p) < 2 (or p < 1/2). See Fig. [5] for the reference 
of the analysis. 
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FIG. 2: The structure of Eq. ©. 
/3 G (0, 1) exists for any p G (0, 1/2). 
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Figure [3] shows numerical results of Nl ■ The parame- 
ters are L = 1, and p = 0.1,0.2,0.3,0.4. The probability 
density g for the fracture points are g(£) = <5(£ — 1/2) 
in (a) and <?(£) = 1 in (b). We limited fracture to 20 
stages at the maximum, and we counted only the inac- 
tive fragments. Each plot is the average of 1000 sam- 
ples. An exact power law Nl is also shown with black 
lines. Power laws fail in larger fragment sizes, as men- 
tioned above. In the panel (b), the cumulative numbers 
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FIG. 3: Numerical result of cumulative number Nl{x) for 
L = 1, and p — 0.1,0.2,0.3, and 0.4, generated by counting 
only the inactive fragments within 20th stage of cascade, and 
averaging 1000 samples each. Each straight line indicates 
the corresponding Nl that follows an exact power law. The 
fracture points are (a) at the middle of the fragments = 
S(x — 1/2), and (b) distributed uniformly = 1. 



also deviate largely from power laws in smaller sizes be- 
cause the number of fracture steps is bounded: some long 
fragments are still active after the last fracture stage in 
the simulation, and they will contribute to raising of the 
number of small fragments if fracture is continued. 

It is noted that there have been many experimental re- 
sults of (3 > 1, but the above model provides only /3 < 1. 



Here we modify the model in order to realize (3 > 1. Re- 
calling the above analysis, we have treated the stopping 
probability as a constant value p. Instead, we set here the 
stopping probability as a function of a fragment size. In 
particular, we give the stopping probability of a fragment 
of size £ as 



(£<X) 
(£>X) 



(5) 



where A is a characteristic length and 7 > is a con- 
stant. It represents an effect that smaller fragments are 
more difficult to experience further fracture. Obviously, 
a fragment becomes inactive whenever its size becomes 
smaller than A, hence the parameter A is the lower bound 
of the fragment sizes. We employ the assumption A <C L 
in the following analysis. 

As above, the cumulative number Nl,\{x), including 
two parameters L and A this time, plays an important 
role in the following analysis. In the same way as Eq. 
©, Nl,\ satisfies the following equation. 

N L ,x(x) 

= M - (J) 7 }/ {^,a(z)+%-«)l,a(x)}<?(£R 

A^ 7 
L 



{N^ tX (x) + N a _ e)LtX (x)} fl (0d& 



(0) 



where we used the approximation A/L ~ 0. (the symbol 
"~" is used only in this sense.) 

A scaling relation N a L^ a x(ax) = N^^xix) is again ob- 
tained. We need another scaling relation for the analysis. 
By the definition of the cumulative number, 



N L ,x(x) = J2 P{£i ■■■ZnL> x){l - px(L)}{l - Px&L)} ■ ■ • {1 -pa (6 • • • t n -iL)} • Pa(6 ■ • • 

n=0 

°° / \ \7 

a:) 7 FT ) kA7 - 

^0 Ui •••£»£/ 

1 



Thus, NL, a \{x) — a J NL t x(x) is derived for x ;§> A and 
a > 0. We guess a power-law form Nl,x(x) — Cx~P, 
and substitute into Eq. ([6]) together with two scaling 
relations, which yields 

f{e M + (i-r 7 }sm = i. 

Jo 

This equation can be solved immediately as (3 = 1 + 7, 
where we note the normalization g(£)d£ = 1. j3 > 1 is 



attained because 7 > 0. A remarkable point is that the 
exponent /3 = 1 + 7 is universal over any probability den- 
sity g governing the fracture points. (Compare with the 
case of a constant stopping probability, where (3 depends 
on g.) 

The coefficient C is A 7 L, derived from the consistency 
of two expressions N L {L) = CL~' 3 = CL~ (1+ ^ and 
Nl(L) = px{L) = (A/L) 7 . Finally, the complete solu- 
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FIG. 4: Numerical results of Nl for L = 1. The parameter values of the stopping probability are A = 10" b and 7 = 0.5, 1, 
and 1.5. Each data set was generated by averaging 1000 samples. Solid lines indicate the corresponding solutions JJJ. The 
probability densities for the fracture points are respectively g(£) = 8{£ — 1/2) in (a), and <?(£) = 1 in (b). 



tion is expressed as 



AV /X\-(l+7) 



(7) 



The calculation is based on x » A; consequently, this 
solution probably breaks down if x < A. 

Numerical results are shown in Fig. 0] where we set 
L = 1 and A = 10~ 6 . Numerically-generated cumulative 
numbers clearly lie on the power-law solution (solid line) 
over a wide range of larger fragment sizes. Also, the 
data points deviate from the power laws in a fragment 
size close to or less than A, as expected theoretically. 

One can straightforwardly extend the model so that 
each fragment breaks into n subfragments at a single frac- 
ture, where n can be either a fixed or random number. 
A fragment size distribution in this case is again like a 
power law; the exponent ft is less than 1 under a con- 
stant stopping probability, and ft = 1 + 7(> 1) under the 
stopping probability as in Eq. (j^J. A special case like 
the Sierpinski fractal is found in Refs. [IH, Qj| without 



pointing out the sensitivity of ft against g. 

Our model claims that a lognormal and power-law dis- 
tribution are similar; the difference is whether the stop- 
ping probability exists or not. Their similarity has been 
supported experimentally. A fragment size distribution 
qualitatively changes according to impact energy 17J (or 
falling height (Isj): it exhibits a lognormal distribution 
under lower energy, and a power-law distribution un- 
der higher energy. These results imply that a lognormal 
distribution and a power-law distribution can possess a 
common origin. Furthermore, the proposed mechanism, 
where a multiplicative stochastic process with random 
stopping produces a power-law distribution, is quite sim- 
ple and general, so it will be applicable to other systems 
than fracture. 
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